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We present the results of the modelling of CO adsorption and catalytic CO oxidation on inhomo- 
geneous Pt(100) surfaces which contain structurally different areas. These areas are formed during 
the CO-induced transition from a reconstructed phase with hexagonal geometry of the overlayer to 
a bulk-like (1 x 1) phase with square atomic arrangement. In the present approach, the surface tran- 
sition is explained in terms of nonequilibrium bistable behavior. The bistable region is characterized 
by a coexistence of the hexagonal and (1 x 1) phases and is terminated in a critical bifurcation point 
which is located at (T c ~ 680 K, p^o ~ 10 Torr). Due to increasing fluctuations, the behavior at high 
temperatures and pressures in the vicinity of this cusp point should be qualitatively different from 
the hysteresis-type behavior which is typically observed in the experiments under ultrahigh vacuum 
conditions. On the inhomogeneous surface, we find a regime of nonuniform oscillations characterized 
by random standing waves of adsorbate concentrations. The resulting spatial deformations of wave 
fronts allow to gain deeper insight into the nature of irregular oscillations on Pt(100) surface. 



PACS numbers: 82.40.Bj, 82.45. Jn,05.45.-a,82.20.Wt 



I. INTRODUCTION 

The processes of surface reconstruction play a central 
role in heterogeneous catalysis A prominent exam- 
ple of a catalyst with the reactivity strongly affected by 
the rearrangement of the substrate atoms is the Pt(100) 
surface. The mechanism of such rearrangements is the 
tensile excess stress due to charge depletion of 5d or- 
bitals of surface Pt atoms 0]. As a result, the bulk-like 
(1 x 1) termination of clean Pt(100) corresponding to 
the square configuration of Pt atoms becomes metastable 
which causes a transformation to the hexagonal (hex) 
surface atomic arrangement. From the point of view of 
the catalytic properties, the crucial importance has the 
fact that such surface reconstructions are strongly af- 
fected by the adsorbates. The CO-induced lifting of the 
hex surface reconstruction and the stabilization of the 
(1 x 1) phase is forced by a gain in the CO adsorption 
energy and proceeds through the nucleation of (1 x 1)- 
islands with high CO coverage 0,0, @. 

The reconstruction of Pt surfaces is a key factor re- 
sponsible for their complex oscillating behavior and pat- 
tern formation Q . This is also supported by the fact that 
the stable Pt(lll) surface where the reconstruction is not 
observed, does not exhibit oscillation properties. Dur- 
ing the reconstruction, the surface of Pt(100) becomes 
inhomogeneous and contains coexisting micrometer-size 
regions of different (hex and square) atomic configura- 
tions. As demonstrated in low energy electron diffrac- 
tion (LEED) studies 0,B@, the propagating waves of 
surface modification can result in highly irregular spatio- 
temporal character of the oscillations. 

Even with a constant oscillation period, the time 
variation of the local LEED intensity of small surface 
spots differs from the temporal behavior of the inten- 
sities integrated over the whole surface area [g]. Such 
a behavior is in contrast to the uniform oscillations 



on Pt(llO) surfaces originating from the global cou- 
ling of various surface patches through the gas phase 
|. It should be noted that the mechanism of the 
global coupling has been extensively studied in a series 
of experimental works and by the theoretical modelling 
[H El El El Q E EDI- In distinction to Pt(llO), 
the irregular oscillations during CO oxidation on Pt(100) 
possibly originate from the propagating chemical reaction 
fronts which can be triggered by defects and other surface 
imperfections @, G3] ■ The kinetic Monte Carlo simula- 
tions performed for the heterogeneous Pt(100) surfaces 
containing mesoscopic hex regions in macroscopic (lxl) 
surroundin g, h ave essentially shown local unsynchronized 
oscillations [18j. These local oscillations disappear on the 
/xm-scale, whereas the experiment demonstrates well de- 
veloped irregular oscillations even on the macroscopic 
mm-scale Q . Another aspect which can lead to the irreg- 
ular temporal character of oscillations is the stong fluc- 
tuations on nm-size catalyst particles with (lOO)-facets 
[19( . As the fluctuations in surface coverages are increas- 
ing with the decrease of the system size [2(| [2l|, they 
cannot satisfactorily explain the irregularities observed 
on single crystal surfaces of Pt(100). 

To gain deeper insight into the nature of the irreg- 
ular behavior, in the present work we consider a mod- 
elling of spatially inhomogeneous Pt(100) surfaces. In 
our approach, we account for the fact that the hex- 
square transformation occurs through a hysteresis and 
has a character of a first-order phase transition. Un- 
der the non-equilibrium conditions, the hysteresis fea- 
tures are described in terms of a bistable behavior which 
is characterized by a coexistence of the hex and square 
phases in a wide range of temperatures and partial pres- 
sures. Such a mechanism allows to obtain in a natural 
way the inhomogeneous surface state, without introduc- 
ing additional separate variables for the description of 
the adsorbate coverage in different phases. Upon the in- 
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crease of temperature and CO pressure, the bistable re- 
gion narrows and terminates at a bifurcation cusp point 
(T c w 680 K, p co w 10 Torr). Due to strong fluctuations 
near this critical point, the behavior at high tempera- 
tures and pressures should be qualitatively different from 
the hysteresis-type behavior observed on Pt(100) under 
ultahigh vacuum conditions 

Furthermore, we address two main questions related to 
the surface inhomogeneities. First, we analyze the role 
of the inhomogeneities in the nucleation and growth of 
adsorbate islands. We find that the difference between 
the CO desorption rates for the hex and square substrate 
geometries is responsible for the trapping of the adsorbed 
CO [|[, a property responsible for the growth of CO is- 
lands during the hex-square transformation. Second, we 
study the oscillation behavior on inhomogeneous surface 
during the reaction of CO oxidation. We show that the 
surface inhomogeneities lead to the Benjamin-Feir insta- 
bility and to a new inhomogeneous oscillating regime. 
This new regime is characterized by the random stand- 
ing waves in spatially separated surface regions. We in- 
vestigate the driving mechanisms of irregular oscillations 
and provide the analysis of the stability of the system 
with respect to nonuniform perturbations. The obtained 
deformations of the chemical wave fronts appear to be a 
direct consequence of structural inhomogeneities typical 
for Pt(100), and can be a possible clue for the explana- 
tion of the irregular oscillating behavior on this type of 
surfaces. 



II. CO-ACTIVATOR MODEL AND ITS 
ANALYSIS 

A fundamental property related to the adsorption of 
CO on low-index Pt(100) surfaces is a transformation of 
the initially stable reconstructed hex Pt surface into the 
bulk-like (1 x 1) structure. The lifting of the surface re- 
construction is activated in the case when the concentra- 
tion (coverage) of adsorbed CO exceeds a critical value 
®CO ~ ^-05 ML. In a wide range of the CO surface cov- 
erages ©co < ©co < ©co (©co ~ °- 3 ML )' tne surface 
contains the coexisting areas of the (1 x 1)- and hex- 
geometry [f|. In the LEED experiments, the range of 
temperatures and CO partial pressures determining the 
phase coexistence is identified by a hysteresis behavior of 
LEED intensities Q. 

It is noteworthy that in the nonequilibrium conditions, 
the hex-(l x 1) coexistence should be described in terms 
of a bistable behavior. To provide a consistent analysis 
of the CO controlled bistability, we introduce a variable 
A = ©co for the description of the CO surface coverage 
and a parameter n which measures a degree of the surface 
reconstruction. The value n = corresponds to a fully 
reconstructed hex surface whereas the maximal value n = 
1 describes the recovered homogeneous (1 x 1) state. 

Within the two-variables approach, the process of the 
CO-induced surface transformation contains two main 



stages. The first stage includes the change of the local 
mesoscopically averaged CO surface coverage during the 
steps of CO adsorption, desorption and CO diffusion on 
Pt surface: 

dA 

— = PA k A (l - (A/A s ) q )s A ~d A -A + D A AA. (1) 

Here p A and k A denote the partial pressure and the im- 
pingement rate of CO, respectively, and A s = 0.5 refers 
to the maximal (corresponding to a saturation) CO cov- 
erage; D A is the CO diffusion coefficient. The factor 
q > 1 models the precursor-type kinetics of the CO ad- 
sorption [22[ I23I]. In our analysis, a typical value q = 3 
for the precursor factor is chosen. 

For the kinetic processes on the inhomogeneous sur- 
face containing a mixture of hex- and (1 x l)-islands, 
we need to consider two different contributions to the 
adsorption/desorption steps in (fT|). These two different 
types originate from the processes CO gas <->COhcx and 
COg a s ^COixi which contribute with the weights (1— n) 
and n to the total rates of surface adsorption and des- 
orption. As a consequence, the final expression for the 
CO sticking and CO desorption coefficients s A and d A 
contain both contributions 

s A = 4(1 -n) + s\n, (2) 
d A = d A (l — n) + d A n, 

where d° A = u exp(- J B°/i?T) and d\ 
Uiexp(—Ej/RT). In the expressions for s A and 
d A , the parameters s A , d A and s\, d\ are the CO 
sticking and desorption coefficients in the hex and 
(lxl) phases, respectively. As is shown in [{|, the low 
desorption rate d\ is a prime factor responsible for high 
CO coverages of the (1 x 1) surface. In contrast, due 
to higher desorption rate on the hex surface, significant 
part of adsorbed CO molecules desorps which leads to 
substantially lower CO coverages in the hex state. 

For the sticking coefficients, we use the values s A — 
0.75 (A = 0.05, hex geometry) and s\ = 0.34 - 0.7 
(A = 0.5 — 0.3, lxl geometry) reported in Ref. 3. In 
the Arrhenius-type form of the desorption coefficients d A 
and d\ we choose the following values for the rate pa- 
rameters: v = 4 x 10 13 s _1 and V\ = 3 X 10 15 s -1 ; 
E° d = 27.5 kcal/mol and E° d = 34.5 kcal/mol 0]. 

Furthermore, in our approach we consider the transfor- 
mation of the substrate state activated by the adsorbed 
CO. In the two- variable model, the time evolution of the 
local mesoscopically averaged surface state parameter n 
is described by the following phenomenological equation: 

— = 7ol . n (l-n)(A/A s ).e(A-A 1 c ) 

-710 • n(l - A/A s ) 2 + nAn. (3) 

Here the first term refers to the lifting of the hex recon- 
struction. The rate of the lifting is proportional to the 
CO coverage a and to the factor 7 in(l — n) which is 
a mean-field form of the effective flux between the hex 
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and (1 x 1)- regions. By introducing the step function 
B(A- Al) (6(£) = 1 for £ > and otherwise), we 
account for the experimental observation that the tran- 
sition to the (lxl) phase does not occur below a crit- 
ical coverage A\ = O^o- ^ n our calculations, we de- 
scribe the transition rate 701 = voi exp(—Eoi/RT) by an 
Arrhenius-type expression where the activation energy 
Eqi ~ 2 kcal/mol corresponds to a small activation bar- 
rier for the transformation from the hex- to the (1 x In- 
state [4], and the prefactor v m k, 10 s _1 gives correct 
time scales for the transformation in the considered tem- 
perature interval [22|, US] . 

The second term in ([3]) refers to the reverse (1 x 1)— >hex 
reconstruction with a rate coefficient 710- Here the fac- 
tor (1 — Aj A s ) 2 accounts for the property that the re- 
construction can occur locally if a small cluster of neigh- 
bouring surface positions (a pair in the simplest case) is 
not occupied by adsorbed CO which is in analogy with 
the assumption considered in [H, [25[ . 

For an inhomogeneous surface, the laplacian term An 
in §5§ originates from the contribution of the interfaces 
between different surface geometries to the total system 
energy [26|, [27j ■ Consequently, the coefficient k describes 
the energy costs of such interfaces a nd is r elated to the 
characteristic interface width L n = y/ ft/701 which is typ- 
ically of the order of few nanometers. We note that the 
characteristic diffusion length of the adsorbed CO is given 
by La = \J Da/joi and is of the order of several microm- 
eters. The fast CO diffusion between the hex- and (1 x 1)- 
areas described by the diffusion term in (TTJ) , controls the 
um size of CO islands which is discussed in details in 
section IHII 

The equations §Q and can be transformed by the 
substitution: 

A = a ■ A s ; x = xLa] cLa = dyt/701; * = <7oi; 
Pa = PAkA/A s ~/ai; r) = k/Da; 7 = 710/701 (4) 

into the following dimensionless form 
da 

TTf = fi{a,n) + A £ a (5) 
at 

—~ = / 2 (a,n) +rjA i: n, 
at 

where 

/1 = p A (l - a)s A -d A -a, (6) 
I2 = n(l — n)a • 0(a — a c ) — 771(1 — a) 2 . 

For a given a ^ 0, the second equation ([5|) leads to the 
two stable steady states: n% = and ri2 = 1 — 7(1 — 
a) 2 1 a. In the limit a — > 1, the second solution approaches 
the value n% = 1 which implies the full lifting of the 
hex geometry by the CO adsorption. In the following 
discussion, we will omit for convenience the tilde signs in 
the notations of dimensionless time and coordinate. 

A remarkable property of the system (0, © is 
a strong dependence of the CO adsorption and desorp- 
tion rates on the local configurational state of the surface 



determined by the parameter n. In analogy to [f|, this 
dependence can can be conveniently analyzed in terms of 
a net sticking probability of CO molecules defined as 

s A = s A - d A a/p A (l - a). (7) 

Specifically, despite a large initial sticking coefficient s ^ 
in the hex regions, the significant contribution of the sec- 
ond dcsorption term in ([7]) results in a suppression of 
which reflects the fact of low CO coverages A < 0.05 ML 
observed on the reconstructed surfaces. In contrast, due 
to the low desorption rates on the (1 x 1) substrate, the 
adsorption of CO will result in large values of s A and 
consequently in high CO coverages in the regions with 
the lifted reconstruction. In the following analysis, we 
will use the net sticking probability as a relevant quan- 
tity for the interpretation of the growth kinetics of CO 
islands on inhomogeneous substrates. 

A. The surface transition hex-(l x 1) on Pt(lOO): 
bistable behavior 

A central feature of the model ([5]), © is the occur- 
rence of a bistable region characterized by the coexis- 
tence of the states m and ri2 with low and high CO 
coverage. In the analysis of ([5|) and ([6]), we choose 
k A = 1.13 ■ 10 5 s _1 Torr _1 ML and perform the calcu- 
lations of the steady states using the rate parameter 
7io = ^10 ex P (—E10/RT) with E\q — 23 kcal/mol and 
U10 = 10 12 s _1 . The large value of E\o is consistent with 
the high activation barrier for the hex reconstruction of 
the clean Pt(100) surface reported in Q. In the Arrhe- 
nius form of the desorption rate parameter d\ of the 
(1 x 1) substrate, we choose the values v\ = lx 10 15 s _1 
and E\ = 37 kcal/mol. These values have been deduced 
in Q from the temperature dependences of the corre- 
sponding LEED intensities which exhibit a hysteresis- 
type behavior. 

Fig. [1] shows the stationary solutions calculated for 
different CO pressures. Here the bistable region where 
the stable hex and (lxl) phases coexist is located in 
the range of pressures between p\ — 3 ■ 10~ 7 Torr and 
p 2 A = 10~ 5 Torr. 

The regions of the stability of different phases together 
with the bistable region are indicated on the phase dia- 
gram (1/T, pco) shown in Fig. [51 We note that the topol- 
ogy of the diagram and the location of the bistable region 
agree well with a LEED diagram obtained by the estima- 
tion of isosteric heats of CO adsorption Q . A remarbable 
feature of our diagram is that the bistable region narrows 
with the increasing pco an d terminates at a critical bifur- 
cation point at p c co w 2 x 10 1 Torr and T c » 680 K. The 
existence of such a cusp point on the diagram is a clear 
manifestation of the critical behavior which is associated 
with a strong increase of near-critical fluctuations in the 
system. We should emphasize that in this high pressure 
range, the direct comparison with the experimental data 
obtained at low CO pressures up to 10 -4 Torr cannot be 
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FIG. 1: Stationary CO coverage a and surface reconstruction 
parameter n versus CO pressure. Here T = 450 K. The un- 
stable solutions correspond to the saddle nodes and indicated 
by dashed curves. 



Pt bulk-like 
(1x1) 




1.5 1.6 1.7 1.8 1.9 2.0 2.1 2.2 2.3 2.4 2.5 
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FIG. 2: Phase diagram (1/T, pco) showing the regions of the 
stability of the hex and (lxl) phases with the indication 
of the bistable region where both these phases coexist on Pt 
surface. 



provided. In view of this, additional experimental stud- 
ies conducted at high pressures and temperatures are re- 
quired in order to verify the existence of the obtained 
critical bifurcation point. For Pt(100), such experiments 
would be of a central importance in bridging a " pressure 
gap" between the hysteresis-type behavior at low pco 
and possible critical fluctuations which should exist at 
high CO pressures. 



III. ISLANDS NUCLEATION ON 
INHOMOGENEOUS PT SURFACE 

In the bistable region, the problem of prime inter- 
est is the behavior of adsorbate molecules on a surface 



which contains micrometer-size areas of different geom- 
etry. To study the nucleation and growth of adsorbate 
islands on such an inhomogeneous substrate, we have to 
consider the full set of the reaction-diffusion equations 
(IH) and ([6]). It should be noted that the strong differ- 
ence between the characteristic nanoscale width of the 
hex/(l x 1) interfaces L n and the micrometer scale of 
the CO diffusion length La results in small values of the 
ratio T] in ((4]) which is of the order 10 -3 -10~ 4 . Further- 
more, in our analysis we choose a typical value r\ = 0.005. 
The equations ([5]) and ([5]) have been solved numerically 
using an implicit two-level finite-difference scheme on a 
discrete finite grid with the Neumann boundary condi- 
tions da/dx — 0, dn/dx = 0. In this absolutely stable 
scheme, the standard finite difference approximations of 
the second-order accuracy for coordinate derivatives and 
the first order accuracy approximation for time deriva- 
tives are applied. 
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FIG. 3: Time evolution of spatial profiles of a and n on in- 
homogeneous Pt surface for pa = 3.3 • 10 -6 Torr, r\ = 0.005, 
7 = 0.5 and T = 450 K. The nucleation and growth of two 
islands with high a-coverage occurs in the lxl surface re- 
gions. The bottom panel shows the behavior of the net CO 
sticking probabilities in the central hex (a; = x c ) and in the 
unreconstructed (x — x~) parts of the substrate. 

To demonstrate the influence of the surface geometry 
on the growth of adsorbate islands, in Fig. [3] we present 
the spatio-temporal evolution of islands with high CO 
coverage on inhomogeneous surface. For simplicity, we 
consider a quasi-one-dimensional case described by the 
spatial coordinate x which corresponds to a stripe-like 
surface patterning. The initial geometry of the substrate 
is characterized by two symmetric regions exhibiting a 
high degree of transformation into the (lxl) phase with 
n = 0.5. These regions of the length 15 /xm are located at 
a distance about 35 fim from the center and surrounded 
by the reconstructed (hex)-areas with n = 0. In the cen- 
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tal region, we initially have a preformed CO island which 
is centered at x = x c and located in the hex-area beyond 
the (1 x l)-patches. Furthermore, in the process of the 
temporal evolution shown in Fig. [31 the CO coverage of 
the central hex-area decreases. This decrease is driven 
by the CO desorption which is dominant due to the high 
desorption coefficient <P A . It is remarkable that the dis- 
appearance of the central CO island is accompanied by 
the simultaneous development of two new symmetric is- 
lands with high CO coverage. These two islands nucleate 
in the (1 x l)-parts of the substrate where the desorption 
is suppressed and CO molecules become trapped due to 
the high CO binding energy on the surface. 

The corresponding difference between the kinetic prop- 
erties of different regions of inhomogeneous substrate can 
be interpreted in terms of the time evolution of the net 
sticking probabilities s A shown in the bottom panel of 
Fig. [31 The predominant desorption of CO from the cen- 
tral hex area is reflected in the initial zero net sticking 
probability s A {x — x c ) which slightly increases with the 
expansion of the (1 x l)-areas on the surface. In contrast 
to this, in the symmetric regions with partially lifted sur- 
face reconstruction, the net sticking probability is high 
due to low desorption rates on these parts of the surface 



in Fig.©. 
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suits in the trapping of CO molecules and in the growth 
of the micrometer-size islands. 

The effect of the CO diffusion on the island growth 
becomes more clear when we perform calculations with 
a small CO diffusion coefficient which corresponds to the 
case shown in Fig. Q] In Fig. [H the decrease of Da by 
the two orders of magnitude corresponds to the addi- 
tional prefactor 10 -2 near the diffusion term in ([S]). The 
comparison with Fig. [3J clearly shows that such a slow 
diffusion leads to a strong localization of the new sym- 
metric CO islands in the (1 x l)-areas and to a significant 
slowing down of the hex-lifting process. 
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FIG. 5: Time evolution of spatial profiles of a and n on in- 
homogeneous Pt surface for pa = 3.3 • 10~ 6 Torr, 7 = 0.5 
and T = 450 K. The CO diffusion constant Da is by one 
order of magnitude smaller then that in Fig. [3] The initial 
spatial distribution of CO is homogeneous with a — 0.05. To 
provide a correct comparison of both spatial profiles, the pa- 
rameter scaling is performed relative to the reference data in 
Fig. [31 The bottom panel shows the behavior of the net CO 
sticking probabilities in the central hex [x = x c ) and in the 
unreconstructed (a; = x_) parts of the substrate. 



FIG. 4: Time evolution of spatial profiles of a and n on inho- 
mogeneous Pt surface for pa = 3.3 • 10 ~ 6 Torr, 7 = 0.5 and 
T — 450 K. The CO diffusion constant Da is by two orders of 
magnitude smaller then that in Fig. [3] To provide a correct 
comparison of both spatial profiles, the parameter scaling is 
performed relative to the reference data in Fig. [3] 

In fact, the obtained simultaneous disappearance and 
nucleation of CO islands in the areas of different geom- 
etry can be characterized as an effective island propaga- 
tion. The basic mechanism responsible for this propaga- 
tion, is the difference between the CO desorption rates 
in the (I x 1) and hex surface parts. Consequently, in 
the (1 x 1) regions where the desorption rate is lower, the 
rapid CO diffusion between the hex and (lxl) areas re- 



To analyze the role of the CO desorption in the is- 
land nucleation, we consider an inhomogeneous surface 
where the initial CO coverage is uniformly distributed 
with a(x) = a,Q = 0.05 (Fig. [5]). This is lower then the 
critical coverage a\ = 0.1 necessary for the hex— ^(1 x 1) 
transition. In the course of temporal evolution, the CO 
coverage begins to increase predominantly in the (1 x 1)- 
regions of substrate which leads to the development of 
two symmetric CO islands. Although the absolute stick- 
ing coefficient in the (1 x l)-areas s A is smaller that the 
hex sticking coefficient s A , the l° w CO desorption rates 
lead to substantially higher net CO sticking probability 
in the (1 x l)-regions, which is illustrated by the time 
evolution of s A in the hex (x = x c ) and (1 x 1) (x = X-) 
regions in Fig. [5] (bottom panel). It should be noted 
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that the obtained strong dependence of the net sticking 
probabilities on the geometry of the substrate is in agree- 
ment with the scenario of CO island growth dynamics 
discussed in In view of this, the difference in the CO 
desorption rates appears to be a central mechanism re- 
sponsible for the CO island nucleation and CO trapping 
on inhomogeneous surfaces. 



IV. MODELLING CATALYTIC CO OXIDATION 

The equations (p} and (J3j) represent in fact a CO- 
activator model where the surface transition hex— >(1 x 1) 
is activated by the CO adsorption. On the Pt(100) 
surfaces, this structural surface transformation is a key 
source for complex spatial self-organization and tempo- 
ral oscillations observed in surface reactions like catalytic 
CO oxidation 0, H| . This is in distinction to the oscillat- 
ing behavior on Pt(llO) surfaces where the oscillations 
occur via a coupling through the gas phase. For Pt(100), 
the basic oscillation mechanism consists in a periodic 
switch between the phases with low- and high reactivity. 
This switch is related to the periodic surface transforma- 
tion between the hex and (lxl) surface structures. We 
note that in contrast to the adsorbed CO which activates 
the hex^(l x 1) transformation in the considered temper- 
ature interval 400-500 K, the coadsorbed oxygen plays a 
role of inhibitor for such a transformation. The inhibiting 
function in this case is based on the fact that the oxy- 
gen predominantly adsorbs in the unreconstructed (lxl) 
phase and reacts with the adsorbed CO. The consequent 
decrease of the CO coverage during the reaction drives 
the reverse (1 x 1)— >hex surface reconstruction. In this 
context, the inhibiting character of the oxygen is a ba- 
sis for the oscillatory behavior in the modelling of the 
surface CO oxidation. 



A. Oscillating behavior 

To study the role of the adsorbed oxygen, we introduce 
the additional equation for the description of the oxygen 
coverage B = Oo on the Pt surface: 

dB , / A B \ 2 

_= PB W„(l-_-_j -r-AB (8) 

where ps = Po 2 an d ks refer to the O2 partial pressure 
and impingement rate, respectively. As the O2 stick- 
ing coefficient in the hex phase is negligibly small, the 
surface state parameter n in the first term restricts the 
oxygen adsorption to the unreconstructed (lxl) parts 
of substrate, which occurs with the sticking coefficient 
s B = 0.3. The mean-field factor (1 - A/A s - B/B s ) 2 
accounts for the condition that two neighbouring surface 
positions not occupied by GO a d s and O a ds a re required 
for the dissociative adsorption of O2. The second term 



in (J8]) describes the decrease of B due to the reaction 
with GO a( is, with the reaction rate r — r exp(—E r /RT). 
The similar reaction term is also included into the kinetic 
equation |T]) for the CO coverage. In our analysis of CO 
oxidation, we choose the values ro = 2- 10 10 s _1 ML _1 and 
E r = 24.1 kcal/mol for the reaction parameters. As the 
diffusion coefficient of adsorbed oxygen is about three- 
four orders of magnitude lower than the CO diffusion 
parameter Da, the adsorbed oxygen is considered as im- 
mobile m. 

After introducing the dimcnsionlcss forms for the oxy- 
gen coverage b = B/B s , oxygen partial pressure Pb = 
Pb^bSb/joiBs and for the reaction rate f — r/701, the 
equation |[8]) can be rewritten in the form 

db 

—= = f 3 (a, n, b) = p B ■ n(l - a - b) 2 - fA s a ■ b. (9) 
at 

In the equation ([5]), the modified function f\(a, n, b) now 
includes the reaction term: fi = PaO- — o,)sa — d,AQ> — 
fB s a ■ b. 

To analyze the stability of the macroscopically homo- 
geneous solutions of ([5]) and ([9]), we performed a lin- 
earization of the system in the vicinity of the stationary 
states £ = (do, n ,bo) T . In this approach, the small de- 
viations 5£ = (5a, Sn, Sb) T from £ can be determined 
from the system of equations 

6£ = G- 5£, (10) 

were G = {gn} is the matrix of the derivatives gn = 

Furthermore, the stability analysis can be reduced to 
the calculation of the Lyapunov exponents Xj (eigen- 
values of G) with the corresponding deviation vectors 
(eigenvectors) <$£■, so that G ■ S£j = (j = 1,2,3). 
From the equations (|10[) . we can easily obtain the expres- 
sions for the time evolution of the eigenvectors 6£j : 5£j = 
c° exp(Xjt) where the constants c° should be found from 
the initial conditions. As a consequence, each arbitrary 
deviation 6£(t) can be represented as 5£ = Y]j c j$£j an d 
in this way is fully determined by Xj . 

In the analysis of the monostable, bistable and oscilla- 
tion states, the Lyapunov exponents have been calculated 
numerically. Fig. [5] shows an example of the oscillating 
behavior in the vicinity of a subcritical Hopf bifurcation 
point. Here, for larger pa, the real part of two complex 
conjugated Lyapunov exponents ReA2 = ReA3 > cor- 
responds to an unstable limit cycle in the phase space 
determined by the vectors £ = (a,n,b) T . Such a cy- 
cle involves several kinetical stages which are based on 
a complex interplay between the CO activator-induced 
hex^(l x 1) transformation, and the inhibitor-caused re- 
verse decrease of the adsorbate coverages resulting in the 
hex surface reconstruction. 

Fig.[6]demonstrates how the periodic oscillations of the 
coverages a and b are related to the surface transforma- 
tion described by the parameter n. The increase of the 
oxygen coverage 6 up to a saturation value 6=1 results 



80 160 
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FIG. 6: Time evolution of dimensionless CO and oxygen cov- 
erages a and b and of the surface state parameter n at various 
CO pressures pco = Pa- Here po 2 = Pb = 5 ■ 1CP 5 Torr, 
k B = 5.6 ■ 10 5 ML-Tbrr _1 s _1 , 7 = 0.5 and T = 450 K. The 
stages (i)-(iii) on the top panel form a single oscillation cycle. 



in a suppression of the oxygen adsorption (stage (i) in 
the top panel in Fig. [5]) . Consequently, the reaction step 
becomes dominating and leads to a decrease of the oxy- 
gen coverage b. During the CO oxidation, the decreasing 
CO coverage a results in the reverse hex reconstruction 
of surface, a process reflected by a slight decrease of the 
state parameter n (stage (ii) in Fig. [5]). The reconstruc- 
tion continues until the CO adsorption begins to prevail 
and the CO-activated hex^(l x 1 ^transformation starts 
again (stage (hi) in Fig. [5J). At this stage, the increase 
of the oxygen adsorption in the (lxl) phase results in 
a repeat of the oscillation cycle. 

As is demonstrated in the bottom panel of Fig. [51 the 
decrease of the CO pressure pa leads to the disappear- 
ance of the oscillations and to a further convergence of 
the system state to a monostable hex phase. 

Fig. [7J presents a bifurcation phase diagram {pa, Pb) 
calculated at T = 450 K with the indicated areas of 
monostable, bistable and oscillating behavior. One can 
see that the oscillating region extends with the increasing 
Pb which supports a key role of the oxygen in the devel- 
opment of the oscillations. Another remarkable feature 
is that the region where the oscillations occur coincides 
with the region of the structural surface transformation 
which is in full agreement with the experimental obser- 



FIG. 7: Bifurcation phase diagram (pco, Po 2 ) calculated at 
T = 450. 

vations by LEED ||. 

On the diagram shown in Fig. the oscillating regime 
disappears in the limit ps — > which corresponds to the 
absence of the adsorbed oxygen on the surface. It should 
be noted that even at extremely low oxygen partial pres- 
sures pb < 10~ 7 Torr, the calculations still give a region 
with low oscillation amplitudes induced by small nonzero 
b. This is in contradiction to the experimental indica- 
tions of the disappearance of oscillations for pb < Pa, & 
property explained in terms of a blocking of the adsorp- 
tion sites by CO Q. Due to a mean-field form of the 
equations used in our modelling, even at low pressures 
p B < 10- 7 Torr we find a small nonzero b which leads to 
the existence of the oscillations. To improve the results 
in this range of pressures, more precise methods have to 
be applied. A possible way to account for the blocking of 
the adsorption sites by CO would be to perform Monte 
Carlo simulations which will require further numerical 
investigations. 

The crucial difference between the oscillations on 
Pt(100) and Pt(110) is related to their spatio-temporal 
behavior. While Pt(110) oscillations arise due to the 
coupling through the gas phase and have a well devel- 
oped regular character, the oscillations on Pt(100) are in 
general irregular. This irregularity is usually explained 
by the inhomogeneous nature of Pt(100) surfaces, where 
the waves of structural transformation develop on struc- 
turally different surface patches at different time inter- 
vals @. Due to a weak spatial coupling via the surface 
diffusion, the resulting integrated oscillation profiles are 
nonuniform and highly irregular Q. 

To gain more insight into the nature of irregularities 
on Pt(100), one needs to study comprehensively the os- 
cillations on surfaces which contain patches of different 
geometries and characterized by different levels of ad- 
sorbate coverages. In fact, such studies can be consid- 
ered in a more wide context, since the modern litho- 
graphic techniques give a possibility to prepare prepat- 
terned surfaces^. Moreover, the control of adsorbate 
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coverages can be achieved by the methods of scanning 
tunneling microscopy and laser induced thermal desorp- 
tion which allow to create artificial micrometer-size mi- 
crostructures on the surfaces |29| . In this way, the char- 
acter of the oscillating behavior can be probed and ex- 
ternally tuned. 



V. CONSEQUENCE OF SURFACE 
INHOMOGENEITIES: CHEMICAL 
TURBULENCE 

To study the effect of inhomogeneities, we consider a 
one-dimensional substrate of a size L x . The substrate 
contains a partially reconstructed central surface region 
C surrounded by the (1 x l)-surface parts where the hex 
phase is fully lifted. For x <E C, the initial surface state 
is described as n(x,t = 0) = n p < 1. On the substrate, 
the region C has the length 2A < L x and is defined as 
\x — x c \ < A where x c is the center of the substrate. In 
our analysis, we choose the following initial conditions 
for the surface coverages: 



b(x,t = 0) = 



bp, 
0. 



xGC 



> A 



; a(x,t = 0) = a p (ll) 



Inside C, the temperature, reaction/diffusion param- 



eters and the initial conditions £ t 



I (2 D , TLj 



respond to the homogeneous oscillating state close to 
the subcritical Hopf bifurcation point. This implies 
that in the case when C will cover the entire substrate 
(2 A = L x ), the time behavior is characterized by the 
homogeneous periodic oscillations of the state parameter 
n(x, t) = n{t) and of the surface coverages a(x, t) = a(t) 
and b(x, t) = b(t) already considered in previous section. 
In contrast, in the case when the oscillating surface region 
C is located inside the unreconstructed surface area, the 
gradients of the adsorbate coverages and of the state pa- 
rameter n near the hex/(l x l)-interfaces at x = x± lead 
to the transition to a highly nonuniform state. To see 
this, in Fig. [5] and Fig. Owe present the evolution of the 
substrate geometry and adsorbate coverages for a case 
when 2A = 40 /xm< L x — 70^m. One can clearly ob- 
serve the development of completely new inhomogeneous 
surface state. In this state, the initial spatial plateau- 
like distribution of b is broken. Moreover, the profiles 
of the oxygen coverage b can be identified as inhomoge- 
neous oscillating standing waves. At the first stages of 
the time evolution (t < 38 sec in Fig. [H]) , these standing 
patterns are localized predominantly within C. Further- 
more, the additional irregular patterns develop near the 
substrate boundaries due to strongly increased b, which 
leads to a final extension of the inhomogeneous oscilla- 
tions well beyond C on the entire substrate (cases t = 110 
and 472 sec in Fig. [8]) . The breaking of the homogeneous 
oscillating regime corresponds to the Benjamin-Feir in- 
stability of homogeneous (fc — > 0) oscillating waves. The 
appearing irregular spatial patterns are characterized by 



a deformation of the wave fronts observed in the ampli- 
tude maps of b in Fig. [TD] These deformations are related 
to a transition to the regime of phase turbulence [30l ]. 
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FIG. 8: Time evolution of spatial profiles of the oxygen 
coverage 6. Here the initial distribution of b is step-like 
with b = bp — 0.63 inside the central region C of a length 
2A = 40fim with x c = 35 /im and 6 = beyond this region. 
Here p A = 6.5 • 10" 6 Torr, p B = 10" 5 Torr, 77 = 0.005, and 
T = 460 K. 
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FIG. 9: Time evolution of spatial profiles of the surface geom- 
etry described by n and of the CO coverage a. The initial n 
distribution is step-like with n — n p = 0.57 inside the central 
region C of a length 2A = 40/xm with x c = 35 \i m and n = 1 
beyond this region. The initial CO coverage a(x) — a p = 0.41 
is homogeneous. Here pa = 6.5 ■ 10 -6 Torr, ps = 10 -5 Torr, 
r) = 0.005, and T = 460 K. 

To deeper understand the origin of the standing waves, 
let us consider the equations (O, © governing the spatio- 
temporal evolution of the system. In the development of 
the wave instability, the property of the CO coverage a to 
adjust the surface state plays a key role. This is demon- 
strated by the temporal evolution of a and n presented 
in Fig. [5] One can clearly see the development of in- 
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homogeneous spatial profiles from the initial uniformly 
distributed a(x) — a p . The decrease of a inside C is pro- 
duced by the initial plateau-like distribution of n in this 
region. As a consequence, the parabolic concave shape of 
a(x) leads to a negative second derivative d 2 a/dx 2 close 
to x = x±. To see the effect of d 2 a/dx 2 < on the wave 
properties, we consider a modified equation for b which 
can be derived in the vicinity of £ p on the basis of ^ 
and ©: 

d 2 b d 2 a 

^2^g^( 2n pPB-rAs)b + f . (12) 

Here /o ~ — 2n p ps(l — a)d 2 a/dx 2 . In the region C, 
the expression for d 2 a/dx 2 can be approximated by a 
parabola: 

d 2 a/dx 2 w - ai {x - x c ) 2 + a , (13) 
with the curvature determined by the parameter a±. 




time (sec) 

FIG. 10: Amplitude map of oxygen coverage b in the (t, 
x) plane in the regime of inhomogeneous oscillations. Here 
p A = 6.5 • 10" 6 Torr, p B = W' 5 Torr, and T = 460 K. The 
left panel shows the development of the central nonuniform 
pattern and the expansion of irregular oscillations near the 
substrate boundaries at the initial stages of temporal evolu- 
tion. 

Furthermore, close to the hex/1 x 1-interfaces, we have 
— ai(x c — x) 2 + a = —a(x) < 0. As a result, in the 
interface region the time-dependent part of the solution 



of (fT2"]) has the form of a standing wave 

b(x,t) = b cos('q(x)t + r) ), (14) 
T)(x) = ^Ja(x)(2n p p B - fA s ). 

The parameters 770 arL d &o can be found from the bound- 
ary conditions b(x±,t) = b > cos(ut). Here b > (i) — 
by cos(ojt) is the oxygen coverage in the (1 x 1) areas 
beyond C which is assumed to be homogeneous for sim- 
plicity. The frequency w is the oscillation frequency of 
the system. With these boundary conditions, close to 
the interface Qx — x c \ ~ A) we have 

b(x,t) — 6> cos(u;t + 5<j)(x,t)), (15) 

where the function 6<p(x,t) ~ ai(A 2 — [x — x c ) 2 )t/A de- 
scribes the deformation of the homogeneous wave profile 
due to a nonzero curvature a\ from (|13p . As the phase 
deformation 6<p(x,t) is controlled by a\, the oscillations 
of a\ (t) during the periodic cycles result in the inhomoge- 
neous wave front oscillating with the constant frequency 
uj which is shown in the right panel of Fig. 1101 

In order to test the Benjamin- Fair instability charac- 
terized by unstable oscillating state with respect to inho- 
mogeneous perturbations, we have calculated the Floquet 
exponents for different wave numbers k. For the analysis 
of the spatially inhomogeneous states, we introduce the 
Fourier transforms = (ak,rik,bk) T and fi{k) 

fii k ) = JT- XI M a ' n ' b ) exp(-ifcr) = f l {a kl n k ,b k ){\&) 

x X 

where the range of the wave numbers k = 2ttI/L x (I 
I N x ) is determined by the periodic boundary con- 
ditions £(x) = £(x + L x ). 

For a state represented by the vector £° = (do, uq, bo) T 
oscillating with a period T and described by and §§§ , 
the analysis of the Floquet exponents can be reduced to 
the solution of the eigenvalue problem with the Jacobian 
Gk = {9i,j(k)} where gij(k) = dfi(k)/d£, k j\^- The con- 
dition X k < implies the stability of £° with respect to 
small inhomogeneous perturbations 5£, k . In contrast to 
this, the opposite condition \ k > corresponds to the 
instability of the homogeneous oscillating state and to 
the appearance of modulation in the system. 

Fig. [TT] presents the behavior of the maximal real part 
of the Floquet exponent X k = maxjReA^,, ReA|, ReA|} 
which occurs under a variation of the length A. One can 
clearly see that the decrease of A leads to a strong in- 
crease of Afe for all values of k (Fig. fTTTa)). The positive 
Afe for nonzero k indicates the formation of the inhomo- 
geneous state presented in Fig. [5] and Fig. [ITJ] More- 
over, Fig. [TTT b) shows that the bulk-like oscillation state 
which develops in an initially homogeneous region C of 
a maximal length A = 35/im, is unstable with respect 
to inhomogeneous perturbations. On such a substrate, 
the uniform oscillations cannot be destroyed by homoge- 
neous perturbations which is reflected in X k (k = 0) < in 
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FIG. 11: (a) Dispersion of the real part of the Floquet expo- 
nent Afc for different values of the half-width A of the recon- 
structed region C. (b) A^ as a function of A for different k. 
Here p A = 6.4 • 1CT 6 Torr, p B = 1CT 5 Torr, and T = 460 K. 

Fig. lllf b). In contrast, due to the positive Afc(fc ^ 0), the 
inhomogeneous perturbations diverge and lead to spa- 
tially modulated oscillatins. 

The principal feature of Afc is its almost dispersionless 
(flat) character which can be clearly seen in Fig. [TIT a). 
Specifically, for all values of A we obtain an initial in- 
crease of Afc at small k, with a further saturation to the 
maximal constant value which occurs for larger k > 0.5 
(Fig. 111( a)). This "flat-band" -character is in strong con- 
trast to the dispersion Afc usually obtained in the mod- 
elling of the global coupling for Pt(110) [10, EH- In the 
latter case, Afc has a distinct well-defined maximum at a 
certain k — ko which corresponds to the appearance of 
spatial clusters with fco-modulation. In our case, due to 
the flatness of Afc at larger k, no distinctly leading mode 
can be formed. Instead, all possible /c-modulations can 
exist in such an instable system. This is the key rea- 
son for a highly irregular form of the standing waves in 
Fig. [5] Although keeping the basic features related to the 
central symmetry, the irregular wave profiles are strongly 
sensitive to the initial state of substrate and to the initial 
spatial distributions of adsorbate coverages, and in this 
sence they are random. In the real systems, the type of ir- 
regular profile can be pinned by surface defects like steps 
and surface imperfections, with possible further random 



changes caused by the noise in the system. 

In fact, the obtained spatio-temporal evolution of 
the irregular oscillations on the substrate is consistent 
with the the development of irregular patterns on spa- 
tially separated patches typically observed on Pt(100). 
It should be noted that the inhomogeneous oscillating 
states and standing waves have been previously discussed 
for CO oxidation on Pt(110) where the main mechanism 
responsible for these_phenomena is the coupling through 
the gas phase [T3, El, E3- In contrast to Pt(110), on 
Pt(100) the inhomogeneities in the surface geometry play 
a central role. Such inhomogeneities lead to the nucle- 
ation of adsorbate islands and to a decisive role of the 
boundaries between geometrically different surface areas 
in the formation of standing waves and irregular oscilla- 
tions. 



VI. CONCLUSIONS 

We have studied the CO adsorption and catalytic CO 
oxidation on inhomogeneous Pt(100) surfaces which con- 
tain structurally different hex and (1 x l)-domains. In our 
model, the structural surface transformation is described 
in terms of a bistability characterized by a coexistence of 
the reconstructed and lxl surface areas. The bistable 
region terminates in the critical cusp point which can be 
approached at higher temperatures and pressures. Due 
to significant increase of the fluctuations, the behavior 
in the vicinity of the cusp point should be qualitatively 
different from the hysteresis-type behavior typically ob- 
served in the experiments under ultrahigh vacuum condi- 
tions. We have also studied the nucleation and growth of 
CO islands on inhomogeneous substrate, a process driven 
by the difference in the desorption properties of CO in the 
hex and lxl phases. We have also performed the anal- 
ysis of the oscillation behavior on the surface during the 
CO oxidation. The surface inhomogeneities lead to the 
Benjamin-Feir instability and to the irregular standing 
waves of adsorbate coverages. The obtained deformations 
of the wave fronts is a direct consequence of structural 
inhomogeneities typical for Pt(100) which allows to gain 
deeper insight into the nature of irregular oscillations on 
this type of surfaces. 
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